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1  Introduction 


In  the  microwave  remote  sensing  of  ocean  surface,  the  use  of  polarimetric  passive 
techniques  has  shown  potential  for  enhancing  the  retrieval  of  wind  speed  and  direc¬ 
tions  [1].  Recent  theoretical  and  experimental  research  activities  have  concentrated 
on  studies  of  polarimetric  thermal  emissions  regarding  the  anisotropic  ocean  sur¬ 
face  assuming  a  smoothly  varying  surface  profile  [1,  2,  3,  4].  However,  under  high 
wind  conditions,  the  presence  of  breaking  water  waves,  foam  patches  and  bubbles 
will  affect  the  polarimetric  brightness  temperatures  of  the  plain  ocean  surface.  The 
significance  of  foam  on  the  ocean  surface  was  recognized  a  long  time  ago  [5],  and 
several  subsequent  experiments  performed  have  verified  its  importance  [6,  7].  Pre¬ 
vious  studies  of  the  foam  contribution  to  the  emissivity  of  ocean  surface  were  based 
on  empirical  formulations  [8,  9]  derived  from  experimental  data.  Although  several 
attempts  at  theoretically  modeling  the  foam  have  been  presented  [10,  11],  it  is  diffi¬ 
cult  to  incorporate  them  with  rough  ocean  surface.  The  more  realistic  modeling  for 
foam-covered  ocean  surface  has  been  proposed  by  Huang  et  al.  [12],  who  consider 
the  sea  foam  to  be  a  layer  with  water  particles  over  a  rough  sea  surface. 

However,  it  is  not  suited  to  model  the  sea  foam  as  the  layer  of  spherical  water 
particles,  since  the  sea  foam  is  dominated  by  water  bubbles  [13].  In  this  chapter, 
we  present  the  theoretical  study  on  the  polarimetric  thermal  emissions  from  foam- 
covered  ocean  surface  based  on  a  composite  volume  and  rough  surface  scattering 
model  using  the  radiative  transfer  theory.  We  model  the  locally  foam-covered  ocean 
surface  as  a  random  layer  with  water  bubbles.  The  small  perturbation  method 
(SPM)  is  used  for  random  rough  ocean  surface,  where  the  bistatic  scattering  is 
calculated  up  to  the  second  order.  The  radiative  transfer  equations  for  foam  layer 
are  solved  using  an  iterative  technique.  The  model  predictions  are  compared  with 
measurement  data  [14]. 

2  Formulations  for  Foam  Emission 

2.1  RT  Equations  for  Foam  Layer 

Sea-foam  is  made  of  spray,  small  water  droplets  and  air  bubbles  which  are  generated 
by  wind  tearing  and  further  processions  such  as  bubble  production,  bubble  down¬ 
ward  entrainment,  and  droplet  produced  by  bubble  bursting  as  illustrated  in  Fig.  (1) 
[13].  Since  the  thermal  emission  from  sea-foam  is  dominated  by  water  bubbles,  we 
simplify  the  sea-foam  as  a  water  bubble  layer  as  shown  in  Fig.  (2).  For  simplicity, 
the  top  surface  of  the  foam  layer  is  considered  to  be  a  flat  surface.  Above  the  foam 
layer  (region  1)  is  a  half  free  space  that  is  labeled  as  region  0  with  e0,  n0.  The  foam 
layer  is  specified  by  a  foam  thickness  di,  the  inner  bubble  radius  R,  bubble  film 
thinkness  5  and  permittivity  ei,  the  fractional  volume  fv  of  bubbles,  the  extinction 
coefficient  Ke,  and  the  temperature  profile  T(z).  The  background  of  the  foam  layer 
is  considered  to  be  free  space  with  e0,/v  The  sea  water  (region  2)  is  in  the  lower 
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Figure  1:  The  generation  of  sea-foam. 
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Figure  2:  The  configuration  of  local  foam  layer  on  wind-driven  rough  ocean  surface. 


half  space  with  permittivity  €2,  salinity  S.  and  physical  temperature  Ta.  The  foam 
coverage  is  denoted  by  F,  thus  the  coverage  of  the  plain  ocean  surface  is  1  —  F. 

In  the  foam  layer,  the  radiative  transfer  equation  [15,  page  229]  is  given  by 

d 

cos  6—I(6,<f>,z)  =  — Re(6,<f))  •  I(9,(fr,z ) 
dz 

+Jdn,7(e,<fi,d,,<j>')-7{6',(j)',z)  +  iT(o,<!>,z),  (1) 

47 T 

where  P(9,  <f>,  6',  </>')  is  the  phase  matrix  which  is  derived  by  using  Mie  theory  in 
Appendix  A.  The  phase  matrix  provides  the  contributions  to  the  specific  intensity 
1(6, 4>,  z)  in  the  direction  (6, 4>)  from  the  direction  (O' ,  $).  Ke  is  the  extinction  tensor. 
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In  this  thesis,  the  specific  intensity  with  the  unit  of  W/m2  is  defined  as 

(|£,l2) 


V 


2Re  (EvEj/) 
L2Im  (EvE*h) 


(2) 


where  rj  =  yj fJ-0/€0  is  the  free  space  impedance.  Tr(0,  <f).  z)  is  the  intensity  of  the 
physical  temperature  and  it  is  written  as  7r(0,  4>  z)  —  Ka(9,  (fi)CT(z).  In  this  expres¬ 
sion,  Ka  is  the  absorption  coefficient  vector,  C  =  Kb/ A2  where  Kb  is  the  Boltzmann 
constant  (I<b  =  1.380658  x  10-23  J/K)  and  A  is  the  electromagnetic  wavelength.  We 
assume  that  the  scatterers  are  bubbles  and  the  absorption  is  isotropic,  thus 

Ka  =  K„[l  1  0  Of.  (3) 

As  mentioned  in  the  previous  section,  the  foam  layer  is  modeled  as  the  composition 
of  spherical  water  bubbles  randomly  distributed  in  the  foam  layer,  therefore  the 
extinction  coefficient  Ke(9,4)  is  a  scalar,  i.e.,  %(9,4)  =  ne.  The  formulations  of 
calculating  /ce  and  na  are  provided  in  Appendix  A.  We  assume  that  the  temperature 
T(z)  in  the  foam  layer  is  independent  of  elevation  and  it  is  equal  to  the  temperature 
of  the  sea  water,  i.e.,  T(z )  =  T0.  Thus 

7T(9,4,z)  =  Ra(9,4)CT0.  (4) 

Define  the  new  specific  intensities  for  0  <  9  <  tt/2  as 

7u(9,4,z)  =  7  (9,4,  z), 

Id(9,4,z)  =  7(n-9,  4,z), 

7 tu(0 -.</>■.  z)  =  7 T(9,4,z), 

1 Td(9,4,z)  =  7t(tt  -  6,  4,  z), 

the  RT  equation  (1)  is  split  into  two: 

cos9^-7u(9,4,z)  =  -k£7  u(9,4,z)  +  [  dQ!=P(9,4,9,,4')-7u(9',4',z) 

dz  J 

upper  2tt 

+  J  dSi'7(e,4’,n-e',<t,')-7i(e',4>',z)+7T«(e,4,), 

lower  2ir 

(9) 


(5) 

(6) 

(7) 

(8) 


-cos  9^-Id(9,4,z)  =  -KeId(9,4,z)  +  [  dQ'  P(tt-9,4,9',4')-Iu(9',4',z) 

dz  J 

upper  27r 

+  J  dn'T(it  -9,4, 7T  -  9',  4')  •  7 d(9',  4',  z)  +  7 Td(9, 4)- 

lower  27t 

(10) 
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We  assume  that  the  thickness  of  the  foam  layer  is  much  larger  than  the  penetration 
depth  of  the  electromagnetic  wave.  Under  this  assumption,  there  is  no  returned 
wave  from  the  bottom  of  the  foam  layer.  Thus  the  boundary  conditions  are: 
on  the  upper  boundary  (z  =  di) 

7d(0,<Mi)=7>,<Mi),  (11) 

on  the  lower  boundary  (z  =  0) 

7U(M,0)  =  0,  (12) 

where  7j(0,  <f>,di)  is  the  specific  intensity  of  the  atmospheric  layer. 

2.2  Solution  of  the  RT  Equation 

To  solve  the  RT  equations  (9)  and  (10).  the  iterative  and  numerical  methods  can  be 
used.  The  iterative  method  is  applied  to  scattering  problems  with  small  absorption 
which  is  the  case  for  the  foam  layer,  while  the  numerical  method  can  be  applied  for 
strong  absorption  problems. 

•  In  the  iterative  method  of  solving  the  RT  equations,  we  consider  the  integral 
terms  as  known  from  the  lower  order  solutions,  thus  the  RT  equations  for  each  step 
are  in  the  form  of  an  ordinary  differential  equation  (ODE).  The  general  form  of  the 
first  order  ODE  is  written  as 

^  +  mv(z)  =  g(z),  (13) 

az 

with  the  boundary  condition  y(z0)  =  y0.  The  solution  of  the  ODE  is 

-  /  f{z')dz' 

e  z°  .  (14) 

The  zeroth  order  RT  solution 

To  solve  the  RT  equations  for  the  zeroth  order,  we  first  ignore  the  scattering  by 
assuming  P  =  0.  Thus  we  obtain 

cosfl-^-T JV.'M  =  (15) 

az 

-cose4-7i0)(«,A2)  =  -/cjf(«.*z)  +  7T.J(^,(t),  (16) 

az 

with  the  boundary  conditions 

7i°y,<M,)=7>,<M1),  (17) 


y(z)  = 


f  f  !iz")dz" 

J  g{z')e:°  dz'  +  y0 


Zo 
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if  (9,4, 0)  =  0.  (18) 

Prom  the  formula  of  the  general  solution  (14)  of  the  first  order  ODE,  the  zeroth 
order  solutions  of  the  RT  equations  (9)  and  (10)  can  be  derived  as  follows: 

-cm  -  1  -  e~zsecdKe 

=  ITu(0,<f>) - ,  (19) 

Ke 


2  _  g-{d\-z)  sec  6 Ke  _ 


Kp 


+  7j(e,<l>,d1)e-(d'-z)seceKe.  (20) 


The  first  order  RT  solution 

Plugging  the  zeroth  order  solutions  (19)  and  (20)  into  the  RT  equations  (9)  and 
(10),  and  defining  the  terms  from  the  zeroth  order  solutions  as  follows: 


7^(0,  <M  =  J  dQ!  P{e,4>,  6'. .<})')■  if  (O',  4>',z) 

upper  27 r 

+  J  d9_'T(d,  0,  7T  -  6',  <t>')  ■  if  (6',  <t>\  z)  +  1tu(0.  <f>), 

lower  2; r 

•  V.  -  (21) 


lrlV:Q,z)  =  J 

upper  27t 

+  j  dQ'T(  7T  -  6,  <j>,  7T  -  6',  4>')  -lf{e',  <t>',  z)  +  lTd(0, 0), 

lower  2tt 


(22) 

we  set  up  the  first  order  RT  equations  as 

COS  0^7,  V,  A,  z)  =  A,  x)  +  7^(9,  A,  z), 

(23) 

-  cosO^j-lj'ie.  a,z)  =  -Kj"\e,A,z)+7?MA,  z), 

(24) 

with  the  boundary  conditions 

7™(9,'Mi)  =  7>,.Mi), 

(25) 

7iI)(9,0,O)  =  O. 

(26) 

Again  by  using  the  general  solution  (14)  of  the  ODE,  the  solutions  of  the 
RT  equations  are  obtained  as 

first-  order 

J  sec 9  <i>.  z')  dz' 

.0 


e 


-Z  sec  0/ve 

) 


(27) 
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z 

7^(0,  <f>,  z)  =  -  sec  0  j 7j](0. 4>„  z^'^^dz'  +  7^(0, 4>,  d1)e 


—disec0Ke  gZsecOKe 


The  specific  intensity  in  the  foam  layer  we  are  interested  in  is  the  up-going  intensity 
at  z  =  d\.  By  integrating  over  the  elevation  in  Eq.  (27),  the  first  order  intensity  of 
the  up-going  wave  at  the  top  of  the  foam  layer  can  be  written  as 


_M,  _  I  _  p-disecBue 

=  iTu{0,4>) - : - 

Ke 


■  sec  0  [  dfl'  P(0,  <£,  0',  <f{)  •  ITu(d',  <p')-2 

J  Kg 

upper  2 tt 

J  _  g-disec&Kc  e~di  sec^Ke  _  ^-d\secOKe\ 

sec  0  sec  0  —  sec  O'  ) 


i  jr/\  1 


+  sec0  J  dQ!  P(6,  <fi,  tt  —  O'.  $)  •  <P') 

lower  2tt 

/j  __  e-disec0Ke  ^  __  e-di{secd+sec$')Ke  \ 

l  sec  0  : . sec  0  +  sec  0'  / 


+  sec  0  /  dft'  P(0,  (j),  tt  -  0',  0')  •  7^(0',  <//,  dj) 


lower  2tt 


1  __  p-di(sec0-fsec0')Ke 

- — - .  (2 

(sec0  +  sec0/)/ce 

Assuming  d\Ke  »  1.  the  Stokes  vector  of  the  thick  foam  layer  is  given  as  follows 


•7T(0.<Mi)  «  iu(0,<f>,d1) 


<!>)-+  [  dQ!  P{9,  (j),  0'.  (j)')  ■  Itu{0'-,  4>')~ 

hup  J  hu 


upper  2tt 


+  J  dtt p{e, <}>,*- e\<j)')-iTd{e',<j>')^  fi- 


lower  27r 


sec  0  +  sec  0' 


lower  2tt 


where  the  first  term  denotes  the  direct  emission  due  to  the  physical  temperature  of 
the  foam  layer,  the  second  and  the  third  term  are  the  scattering  of  emission  by  the 
scatterers  in  the  foam  layer,  and  the  fourth  term  is  the  scattering  of  the  atmospheric 
emission.  Notice  that  the  up-going  and  down-going  specific  intensities  due  to  the 
physical  temperature  in  the  foam  layer  are  direction  independent  and  they  are  given 
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(31) 


by 

'1 

Itu(8,  <l>)  =  -7r</(0)  <f>)  =  KaCT0  -  naCT0  ^ 

.0 

2.3  Foam  Coverage 

Let  the  foam  coverage  be  F,  thus  the  total  brightness  temperature  is  [8] 

TTu  (0,  <J>,d1)  =  F-  Tf  (0,  <t>,  0)  +  (1  -  F)  ■  rf  (0, 0, 0),  (32) 

where  T^F(0, 4 >,  c?i)  is  the  emission  of  the  100%  foam  as  in  Eq.  (30),  and  T^S  (0, 4>,  0) 
is  the  emission  of  the  plain  ocean  surface  with  the  consideration  of  the  reflection  of 
atmospheric  emission  that  will  be  discussed  in  the  following  sections.  Notice  that, 
at  this  stage,  we  only  consider  the  emission  at  the  zero  elevation  height  [z  =  0). 
We  need  to  consider  atmospheric  attenuation  and  radiation  if  we  calculate  the  total 
brightness  temperature  at  the  height  of  the  radiometer  (z  =  d2).  The  empirical 
formula  of  the  foam  coverage  F  as  the  function  of  the  wind  speed  Uw.  the  polar 
angle  0,  and  the  operating  frequency  /  was  provided  by  Stogryn  [8]  and  later  used 
by  Pandey  [9]  among  others.  The  foam  coverage  is  expressed  as  [8] 

F  =  b0  +  biUw  +  b2Ul,  (33) 

where  the  coefficients  b0,  b\  and  b2  are  frequency  dependent,  and  they  are  given  by 

£>o  =  1.707  x  10-2  +  8.560  x  10-4/  +  1.120  x  10-5/2, 

bx  =  -1.501  x  10-2  +  1.821  x  10'3/  -  4.634  x  10“5/2, 

b2  =  2.442  x  10“4  -  2.282  x  10 "6/  +  4.194  x  10 ~7f.  (34) 

In  (33)  and  (34),  the  units  of  the  wind  speed  and  the  frequency  are  m/s  and  GHz, 
respectively. 

3  Thermal  Emission  from  Plain  Ocean  Surface 

In  local  regions  without  foam,  the  thermal  emission  from  the  ocean  surface  is  the 
sum  of  the  reflection  of  the  atmospheric  emission  and  the  thermal  emission  from  the 
plain  ocean  surface,  i.e., 

7f(«,<M)=  /  0)+7f(«>,0),  (35) 

lower  27 r 

where  R  is  the  reflection  matrix  of  the  rough  sea  surface,  Ij  (0',  <p\  0)  is  the  thermal 
emission  of  the  atmosphere,  and  Iu  is  the  thermal  emission  of  the  plain  ocean  surface. 
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The  details  of  calculating  the  atmospheric  emission  1^(6',  fa,  0)  will  be  discussed  in 
Section  3.2. 

Considering  that  the  ocean  water  is  in  thermal  equilibrium,  the  Stokes  vector  is 
related  to  the  emissivity  e  by 

=  (36) 

where  Ts  is  the  physical  temperature  of  the  rough  ocean  surface.  By  Kirchholf ’s  law, 
the  emissivity  vector  e(0, fa)  is  the  intensity  of  incident  wave  with  unit  amplitude 
minus  the  total  intensity  of  the  reflected  waves, 


1(6,  <t>)  =  It  -  f 1,(6, 4>]  e„9,:  0)dCl,, 


(37) 


.where 


(1 Evi\ 2)  1 

'  1  ‘ 

1 

(1^4 

1 

V 

2R  e(EviE*hi) 

0 

.  21m  {EviE*hi)  . 

.  0  . 

(38) 


The  amplitude  of  the  v  and  fi-polarized  incident  electric  fields  is  unity,  and  /r  is  the 
reflection  Stokes  vector  of  the  plane  wave  with  incident  angles  6t  and  fa. 


3.1  Stokes  Vector  of  Reflected  Wave 

We  apply  the  small  perturbation  method  (SPM)  to  calculate  the  reflection  matrix  R 
of  the  rough  surface  and  then  calculate  the  reflection  Stokes  vector  Ir.  The  reflection 
Stokes  vector  can  be  written  as 

7r(9,fa0i,fa;0)  =  Ti(e,fa9i,fa)-7i{9i,fa:0),  (39) 

where  R  is  the  reflection  matrix  for  the  Stokes  vector. 

In  the  zeroth  order  SPM  solution,  the  scattered  field  is  specular  and  is  equivalent 
to  the  flat-surface  scattering  problem;  thus  the  reflection  coefficients  of  the  zeroth 
order  solution  are  the  Fresnel  reflection  coefficients.  Using  the  ensemble  average,  the 
second  order  solution  from  SPM  is  also  specular.  Therefore  we  call  both  the  zeroth 
and  second-order  terms  coherent.  The  averaged  field  of  the  first  order  solution  by 
SPM  is  zero,  hence  it  is  incoherent.  By  including  solutions  up  to  the  second  order 
using  SPM,  the  reflection  matrix  is  the  sum  of  coherent  and  incoherent  parts,  R 
and  /?,  respectively,  so  that 

I  =  lC  +  lt.  (40) 

Since  the  scattered  field  from  the  zeroth  and  second  order  SPM  solutions  is 
specular,  the  coherent  reflection  matrix  of  the  Stokes  vector  can  be  written  as 

7l  (i 9 ,  fa  7r  -  0U  fa)  =  [R-j]  <5 (cos  9  -  cos6i)S(<fi  -  fa) 
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i?ll 

^2 

R\ 4 

rc21 

i?22 

■^23 

i?24 

•^32 

^33 

■^34 

^2 

R%z 

^4 

S(cos6  —  cos  9i)6((f>  —  <pi). 


Note  that  the  subscripts  i  and  j  are  associated  with  the  scattered  and  incident 
components  of  the  Stokes  vector,  respectively.  The  element  ( i,j  =  1,2, 3,4)  is 
related  to  the  reflection  coefficient  Rap  with  subscripts  a,  ,6  =  v,  h.  where 

&vt s  _  Rvv(@i  Rvhi@ii<Pi)  &vi  /  a  n \ 

.  Ei,  J  ~  A)  /&(»„*,)  J  [  Ehl  J  ’  W 

and  the  subscript  v  and  h  represent  vertically  and  horizontally  polarized  waves, 
respectively.  By  writing 

r  f 

T  =  i  1^,1 

“  n  2R  e{E‘„Et.} 

.  2Im{£«,JSg} 

and  expressing  E%  and  Efls  in  terms  of  E^  and  Efd  using  Eq.  (42),  it  can  be  shown 
that  _ 

'A  =  [^]  '  7i,  (43) 


where 


Ay]  = 

|Rf,|2  |fl;.„|2  MRX)  — Im  (E%RH) 

|fl£„l2  |fl«J2  MRi,Ru%)  -im  (flJ.flS,) 

!MW)  2Re(fif„BS)  Re(flf,.fl;s  +  fl;,,fl;1a  -Im  (fl(,.fl?,5,  -  fl;:,,fl?,y 
21m («;.„/?£.)  2Im(^„RS)  Im(RJ.„fiK  +fl5„fl£)  Re  -  R^.)  _ 

(44) 

and  the  incident  Stokes  vector  /,  is  defined  as  in  Eq.  (38).  Using  SPM,  (a,  (3  = 
u.  ft)  can  be  obtained  as 


I  pc  |2 


Re(i?^i?^) 

Re(i?L^) 


R‘0V>  +  fl'2,’  fl^  +  fl'? 

p(0)  _L  E>(°)  4.  ft(2) 

T  -H/,,,  *lhh.  *  rihh. 


where  the  zeroth  order  R^  is  the  Fresnel  reflection  coefficient  of  a  flat  surface,  which 


is  given  by 


p(0)  _  e2^1  z  ~  Clkiz 

vv  -  pklz  +  eyk2z' 

n( 0)  _  fcl z  ~  frzz 

hh  klz  +  &2 2  ’ 

p(°)  _  p(°)  _ n 

^vh  “  “  U* 
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The  second  order  reflection  coefficient  R^p  is  given  by 

oc  oo 

R%  =  I  I  dkxdkyW(kxi-kx,kyi-  ky)fi%  (47) 

—  OO  -OO 

where  W(kxi—kx,kyi—ky)  is  the  spectral  density  function  of  the  rough  ocean  surface, 
and  is  the  second-order  scattering  coefficient  as  in  [1]. 

The  incoherent  reflection  matrix  R  (9,  <6;  61,,  0j)  can  be  calculated  by  considering 
the  first  order  SPM  solution.  By  expressing  the  scattering  fields  as 

Ev$  _  fvv{Qi4,\®ii4>i)  fvh{9,6'i9i,6i)  Evi  ('48') 

Ehs  J“  r  [fhvfrbtuh)  fhh{0,<t>\0i:<t>i)\[Ehi\ ’  1  ' 

where  f 0,3(9. 6\  9X,  6i)  is  the  polarimetric  scattering  coefficient  (a.  (3  =  h.  v),  the 
reflection  matrix  can  be  calculated  by  [16] 


where 


L  (9.6:  &i,  <?i) 
[  (\fw\2) 

(IM2) 


R-A^fL’ 


<l/„|2}  (\U\‘)  Re  ((/„./;„))  -Im  <(/„„«,)) 

(|A„I2)  (IA*I2)  Re  ((/,„/,;„))  -Im  ((A,/,;,,)) 

2R e  ((/„/;..))  2R e  {(/,.„/;„))  Re  -Im  -  /,.*/»%)> 

21m  ((/„/;»))  21m  Im  <(/„/4  +  A*/to)>  Re  ((/,./»  -  A*/;,)) 


and  ^4  is  the  illuminated  area.  The  ensemble  averaged  product  of  the  scattering 
coefficients  is  related  to  the  polarimetric  bistatic  scattering  coefficient  as  following: 

i  a  4n  ft,  4>i)f’m(9,  M,<M) 

Wft  *  «  = - 7^: - •  <51) 

From  expression  (49),  (50)  and  (51),  it  can  be  easily  derived  that 


^ Yvvvv 

lyhvh 

Rnlwh 

“Im7 ywh 

cos  9i 

^fhvhv 

^hhhh 

lm')hvhh 

4tt  cos  9 

2Re%vhv 

-  21m 'Yvvhv 

2Re7 ihhh 
^lmthhh 

ii®  ( 'Yvvhh  d"  ')vhhv) 

im  (rivhh  + 

Im  ( ')VVfih  7t>  hhv) 

)vvhh  ~  ^fvhhv) 

li^u(9,6-,0u6i)  = 


A-Kk\  cos2  9rQ0tlu(9, 6\  6i,  6iW(9, 6]  9U  fa ) 
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where  Tapnu  is  the  coefficient  with  explicit  form  given  in  [1].  Finally,  by  knowing 
the  coherent  reflection  matrix  R  from  Eq.  (41)  and  the  incoherent  reflection  matrix 

R  from  Eq.  (52),  we  obtain  the  total  reflection  R  =  R  4-  R  [Eq.  (40)].  Therefore 
the  reflection  Stokes  vector  IT  is  obtained  as 

7r(MAA;0)-(lf  +  S)*7,,  (54) 

and  the  emissivity  can  be  calculated  from  Eq.  (37)  as 

e(6,  <j>)  =  7  i-  J  (Tl  +T?y  lidQi 

•n/2  2ir 

=  Ii  —  J  sin  OidOi  J  d<f>i 

o  o 

3.2  Reflection  of  Atmospheric  Thermal  Emission 

We  write  the  term  representing  the  reflection  of  the  atmospheric  thermal  emission 
from  the  rough  surface  as 

7r(M;0)  =  J  dn,=R{d,(P,7r  (56) 

From  Eq.  (40)  we  know  that  the  reflection  of  the  rough  surface  can  be  written  as 
the  sum  of  the  coherent  and  incoherent  parts,  therefore 

d>;  0)  =  J  dQ!  (IC  + 1')  -7rfV,  4> 0).  (57) 

Thus  the  total  Stokes  vector  in  Eq.  (35)  can  be  written  as 

7f(M;0)  =  J  dQ.'  (RC  +  Td^lAd(e',o'-Qi) 

+^TS  hi  -  J  dQt  (Ti  +  7i)  ■  7*1  .  (58) 

Noticing  that  dO!  =  dO.,  for  the  down-going  Stokes  vector  from  the  atmospheric 
layer,  therefore  it  can  be  shown  that 

7f  (0,  $>;  0)  =  —-T.h  -  T,  J  dfi,  (ff  +  f  )  •  7(  -  &\0)  ■  (59) 

Converting  the  Stokes  vector  to  the  brightness  temperature,  the  brightness  temper¬ 
ature  of  the  plain  ocean  surface  can  be  written  as 

o)  =  -^7f(«,*0) 

=  T,  -  J d(l,  (jP+tf).  [T,  -  A; 0)] ,  (60) 
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—  -  —A  A2  -A 

where  =  TsIi,  and  Td  (0j,<&;O)  =  — /d  (6>t-,  <pu  0)  is  the  down-going  brightness 

Kb 

temperature  of  the  atmosphere  at  the  ocean  surface.  The  down-going  specific  inten¬ 
sity  ld(0j,  ;  0)  from  the  atmosphere  can  be  found  from  Eq.  (77)  and  its  approxi¬ 
mation  can  be  found  in  Eq.  (89)  as 


d2  -sec  efna(z")dz" 

1 d{0,<j>,di)  =  sec  9  I  Ra(z')CT(z')e  dz' 

di 


CT(z0)  —  1-e  -i 
k  n 


—  seed  J  Ka(z")dz" 


where  T(z')  is  the  temperature  profile  of  the  atmosphere,  and  Ka  is  the  absorption 
coefficient.  The  median  elevation  z0  can  be  calculated  by  solving  the  following 
equation: 

j  KaWW -If  K.tfdz'.  .  (62) 

. .  ;  v\  . .  .  ..  d\ :  ..  *  vt;  :: 


3.3  Power  Spectrum  of  Rough  Ocean  Surface 

The  ocean  surface  spectrum  applied  in  this  thesis  was  proposed  by  Durden  and 
Vesecky  [17].  This  surface  spectrum  is  based  on  experimental  data  fitting  and  thus 
it  is  an  empirical  model.  The  Durden- Vesecky  surface  spectrum  is  given  by 

(63) 

where 

<3>(fc,  4>)  =  (l  +  c(l  -  e_sfc2)cos2^  .  (64) 


if  0  <  k  <  2, 
if  k  >  2. 

The  wind  friction  velocity  u*  can  be  found  from  the  following  equation 
it,  h 

Uh  =  04  °S  6.84  x  10_5/u*  +  4.28  x  10~3ti2  _  4.43  x  lO"4’ 


-a74(u^)2 

/  bku*  \al°gio(fc/2) 
Vs+7/:2) 


(65) 


(66) 


where  Uh  is  the  wind  speed  in  the  unit  of  m/s  at  the  elevation  height  h  in  meters 
above  the  mean  ocean  surface,  k  is  the  ocean  surface  spatial  wavenumber,  <j)  is  the 
azimuthal  angle  with  respect  to  wind  direction,  g0)  o,  6,  g.  7,  and  s  are  constants 
with  the  values  of  ao  =  0.008,  a  =  0.225,  b  =  1.25,  g  =  9.81,  7  =  7.25  x  10-5,  s  = 
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1.5  x  10-4.  The  parameter  c,  which  serves  as  the  coefficient  for  the  azimuthal- 
dependent  term  in  the  spectrum,  is  given  by 


where 


2(1  —  R)  /  (1  +  R) 

C~  1  —  D 

(67) 

n  0.003  +  1.92  x  10-3t/i2.5 

jfx .  n  « 

(68) 

3.16  x  10-3C/i2.5 

Jdkk7  S(k)  e-(fc/89-44)2 

D  =  0 

^  OO 

(69) 

/  dkk2  S(k) 

o 


For  example,  for  the  wind  speed  Uh  =  12  m/s  at  height  h  =  19.5  m,  it  can  be 
calculated  that  it*  =  0.46388  and  c  =  0.65139.  To  consider  certain  hydrodynamic 
effects  of  the  ocean  waves,  we  multiply  the  parameter  c  in  the  spectrum  density 
function  W(k,4> )  by  (1  -  d0  cos^),  where  the  parameter  d0  is  determined  by  data 
matching.  The  hydrodynamic  modulation  was  also  modeled  differently  by  multiply¬ 
ing  the  ocean  surface  spectrum  with  a  parameter  h!  based  on  the  slope  of  the  long 
waves’ [17,  18].’  The  modulated  spectrum  is  written  as' 


W(k,<j>,Sx)  =  tiW(k,<i>),  (70) 

where  Sx  is  the  slope  of  the  large-scale  waves  on  the  ocean  surface  and  h'  is  calculated 
as  , 

,/ _  /  1  -0.5sgn(5x)  if  \SX/SU\  >  1.25,  frn 

~  \  1-0.4  Sx/Su  if  \SX/SU\  <  1.25,  1  j 

where  Su  is  the  rms  upwind  surface  slope  which  can  be  calculated  by  using  Eq.  (110). 
The  hydrodynamic  modulation  of  the  ocean  surface  spectrum  using  the  parameter 
h'  is  useful  for  the  two-scale  model  of  the  ocean  surface.  A  detailed  study  can  be 
found  in  Section  (5). 


4  Radiative  Transfer  Equations  for  Atmosphere 

The  atmospheric  contribution  to  the  brightness  temperature  of  the  ocean  surface 
must  be  taken  into  account  since  (1)  the  airborne  radiometer  is  usually  at  a  very 
high  altitude,  hence  the  accumulated  thermal  emission  along  the  path  from  the 
ocean  surface  to  the  radiometer  may  be  significant,  and/or  the  attenuation  for  the 
brightness  temperature  propagating  from  the  ocean  surface  up  to  the  radiometer 
cannot  be  negligible,  and  (2)  there  may  be  a  significant  amount  of  down-going 
thermal  emission  from  the  atmosphere  being  reflected  by  the  ocean  surface.  In  clear 
air  conditions,  the  main  concerns  about  the  atmosphere  for  passive  remote  sensing 
are  the  atmospheric  emission  and  attenuation  due  to  the  contributions  from  gaseous 
oxygen  (O?),  water  vapor  (H2O),  and  suspended  water  droplets  (hydrosols)  [19].  In 
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Figure  3:  The  atmospheric  layer  above  a  foam-covered  ocean  surface. 

the  adverse  conditions,  cloud  and  rainfall  need  to  be  addressed  [19.  20].  In  clear  air 
and  at  microwave  frequencies,  the  electromagnetic  wave  scattering  by  atmospheric 
•  gases  can  be  ignored  [21],  thus  the  radiative  transfer  (RT)  equations  that  are  used 
to  model  the  wave  propagation  in  the  atmosphere  reduce  to  uncoupled  first-order 
differential  equations. 

The  atmosphere  can  be  modeled  as  an  inhomogeneous  layer  with  the  extinction 
coefficient  ne(z),  the  absorption  coefficient  Ka(z),  and  the  temperature  profile  T(z) 
in  terms  of  the  height  2  as  shown  in  Fig.  (3).  We  assume  that  the  scatterers  such  as 
water  vapor,  droplets  and  gaseous  oxygen  are  small  thus  the  scattering  is  ignored 
(h~e  =  ns  +  Ka  Ka).  Therefore  the  radiative  transfer  equations  for  the  specific 
intensity  Iu  and  Td  have  the  following  simple  form  [15,  22] 

cos6ylAu(e,c i>.z)  =  -Ka{z)lAu{.e,4>,z)  +  RaCT{z\  (72) 

r\ 

-COS0^7>,<M  =  -Ka(z)7*(d,<f>,z)  +  RaCT(z),  (73) 

where  C  =  AT# /A2,  Kb  is  Boltzmann  constant  and  A  is  the  electromagnetic  wave¬ 


length. 

Assuming  that  there  is  no  thermal  emission  from  the  upper  space  at  2  =  d2l  the 
boundary  conditions  for  Eqs.  (72)  and  (73)  are  given  by 

7?(S,'Mi)  =  7l(0,<Mi),  (74) 

7iV.  (75) 

where  Iu(9,  <f>,di)  denotes  the  emission  from  the  lower  boundary  at  z  —  d\  (the 


foam-covered  ocean  surface).  From  Eq.  (14),  the  solution  of  the  RT  equations  (72) 


Radiometer 


Atmosphere 


Foam-covered  ocean  surface 


->-  x 
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and  (73)  are  obtained  as 


d2 

seed  J  Ra(z')CT(z') 

di 


d2 

-s ec6  f  Ka{z")dz" 

e  *'  dz' 


d2 

-secO  f  Ka(zf)dzf 

+7l(e,<f>,d1)e 

z' 

<*?  -sec 6  f  Ka{z")dz" 

7j(0, 4>,  di)  =  sec ej  Ra(z')CT(z')  e  d>  dz! . 

di 


(76) 

(77) 


The  specific  intensity  of  the  down-going  wave  in  Eq.  (77)  can  be  calculated  for  any 
polar  angle  9  from  the  absorption  and  temperature  profiles.  From  the  equation  of 
the  down-going  wave  in  Eq.  (77)  we  find  that  the  assumption  that  down-going  wave 
at  height  d2  is  zero  as  in  Eq.  (75)  is  due  to  the  fact  that  the  absorption  coefficient  Ka 
at  z  =  d2  is  very  small.  In  the  next  sections,  we  will  see  from  the  numerical  results 
for  the  US  Standard  Atmosphere  [23]  that  the  value  of  the  absorption  coefficient 
is  negligible  for  altitude  larger  than  10  km.  In  the  RT  theory,  if  the  scattering  is 
ignored  for  the'  propagating  wave,  the  absorption  coefficient  ka  is  two  times  the 
imaginary  part  of  the  complex  wavenumber, 

Ka  —  21m  {k}  =  21m  {  (lO~6Ar  +  l)  fc0}  ,  (78) 


where  k0  is  the  wavenumber  in  free  space,  and  N  is  the  complex  refract ivity. 
In  the  following  sections,  we  will  use  Liebe’s  millimeter-wave  propagation  model 
(MPM)  [24]  to  calculate  the  complex  refractivity  in  the  atmosphere. 

.  The  extinction  coefficient  Ke  is  the  sum  of  the  scattering  coefficient  ks  and  the 
absorption  coefficient  of  the  background  kq.  Quantitatively,  we  find  that  ks  Ka, 
therefore  we  use  na  to  approximate  Ke.  The  dominant  scattering  in  the  atmosphere 
is  due  to  water  vapor  and  suspended  water  droplets  (hydrosols).  The  scattering  and 
absorption  coefficient  are  given  by  [15,  page  157-158] 


ks  =  2 fvk4a3  \y\ 2 , 

,  _  f  I  3e  '2 

JvK- 


es  +  2e 


(79) 

(80) 


where  fv  is  the  fractional  volume  occupied  by  the  water  particles,  a  is  the  radius  of 
a  water  particle,  and 

y  =  - (si) 


es  +  2e 


For  water  droplets  in  the  atmosphere,  the  typical  values  are  a  ~  5  /im,  fv  ~  5  x  10-7, 
and  es  ~  30  +  z40  for  /  =  20  GHz.  Therefore  ks  ~  10-12m_1  and  na  ~  10-5m_1, 
thus  ks  <C  Ka  and  we  can  ignore  ks  and  let  Ke  &  na.  The  gaseous  oxygen  (O2), 
water  vapor  (H2O),  and  suspended  water  droplets  in  the  atmosphere  are  considered 
as  the  principal  absorbers  in  moist  air  [19]. 
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4.1  The  Millimeter- Wave  Propagation  Model 

The  millimeter-wave  propagation  model  (MPM)  model  developed  by  Liebe  [24]  can 
be  used  to  calculate  the  complex  refractivity  N  that  is  related  to  the  complex 
refractive  index  n  by  N  =  106(n- 1).  By  writing  the  refractive  index  as  n  =  n'+in", 
the  effective  permittivity  of  moist  air  is  found  as  following: 

ee  =  (n'  -P  in")2  =  (n,s  —  n"2)  +  i  (2n'n") ,  (82) 

where  k0  is  the  wavenumber  in  free  space.  Thus  the  complex  wavenumber  is  given 

by 

k  =  LO\J ft,eGe  “ —  k0Tl) 

assuming  (ie  =  fi0.  The  measurable  parameters  of  atmosphere  are  (1)  barometric 
pressure  P  in  kilopascal  (lkPa  =  lOmbar),  (2)  temperature  T  in  degrees  Kelvin 
(K),  (3)  relative  humidity  RH ,  and  (4)  mass  concentration  in  grams  per  cubic  meter. 
The  four  measurable  parameters  can  be  converted  to  the  internal  variables  which 
are  useful  in  the  MPM  model  as  shown  in  Appendix  B. 

The  typical  barometric  pressure  profile  is  given  in  [23]  and  it  is  plotted  in  Fig.  (4). 
As  an  example,  we  consider  the  frequency  /  =  19.35  GHz  and  use  the  US  Standard 
Atmosphere  1976  temperature,  barometric  pressure  and  humidity  profiles  to  cal¬ 
culate  the  complex  wavenumber  k  for  ah  electromagnetic  wave  propagating  in  the 
atmosphere.  The  numerical  result  of  the  absorption  coefficient  is  plotted  in  Fig.  (5), 
where  the  real  and  imaginary  parts  are  plotted  separately  in  terms  of  the  wavenum¬ 
ber  in  free  space  k0. 

4.2  Attenuation  and  Emission  of  Atmosphere 

Once  we  obtained  the  key  parameter  —  the  absorption  coefficient  na  from  the  MPM 
model,  the  thermal  emission  of  the  standard  atmosphere  and  its  attenuation  for  the 
propagating  waves  can  be  calculated  by  using  the  RT  solutions  (76)  and  (77).  The 
first  term  in  Eq.  (76)  represents  the  up-going  thermal  emission  of  the  atmosphere 
and  can  be  written  as 

-see#  f  Ka(z")dz " 

7 *  (0,  </>,  d2 )  =  sec 6  J  R,a(z')CT(z')e  ■«'  dz' .  (84) 

dl 

The  down-going  wave  in  Eq.  (77)  represents  the  down-going  thermal  emission  of  the 
atmosphere, 

z' 

djl  —seed  J  Ka(z")dz" 

7i(e,<t>,d1)  =  sec6jRa(z')CT(z,)e  *  dz'.  (85) 

rfi 

The  up  and  down-going  thermal  emissions  can  be  calculated  using  the  result  of 
absorption  and  temperature  profiles  from  the  MPM.  We  can  show  numerically  that 
they  are  very  close  to  each  other  quantitatively. 
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Figure  4:  Temperature  and  barometric  pressure  profiles  in  the  US  Standard  Atmo¬ 
sphere  1976. 
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Real  Part  of  k-1 


Imaginary  Part  of  k 


Figure  5:  The  complex  wavenumber  of  electromagnetic  wave  with  /  =  19.35  GHz 
propagating  in  the  atmosphere  in  terms  of  the  wavenumber  in  free  space. 


The  up-going  thermal  emission  as  in  Eq.  (84)  can  be  written  as 

,  d2 
p  -sec 0  J  Ka(zn)dz " 

7ui(0,<M2)  =  Cr(*o)sec0  J  Ra{z')e  *'  dzr,  (86) 

d\ 

where  is  between  d\  and  d2-  The  down-going  wave  as  in  Eq.  (85)  can  also  be 
written  as 

z' 

^2  -sec0  J  Ka{z”)dzn 

7>,<Mi)  =  Cr(<)sec0_/ *.(.!:>  «  dz',  (87) 


where  z'Q  is  between  d\  and  d2  but  different  from  z0.  It  can  be  shown  that 

<f2  s' 

d?  -sec 0  f  Ka(z")dz"  d2  -sec 6  J  Ka(z")dz" 

/  Ka{z')e  z'  dz'  =  /  Ka{z')e  dl  dz1 

d  i  d\ 

fd2 

-seed  f  Ka{z,f)dzft 

1  -  e  dl 


Since  the  temperature  profile  T(z)  does  not  change  much  for  the  entire  integral 
path.  i.e.  T(z0)  «  F(z'), 


lh(0,  4>,d2)  «  />,  <Mi)  «  CT(z0)  -  1  —  e 

Kn 


—  sec 0  J  Ka(z”)dzfi 


where  z0  is  the  median  elevation  as  defined  in  Eq.  (62).  Figure  (6)  shows  the  numer¬ 
ical  result  of  the  down-going,  up-going  and  approximated  brightness  temperatures 
by  carrying  out  the  integration  for  Eqs.(84)  and  (85)  numerically  or  using  approxi¬ 
mation  formula  Eq.  (89)  for  the  standard  atmosphere  at  /  =  19.35  GHz,  d2  =  30  km 
and  di  =  0. 

In  the  second  term  of  the  up-going  wave  in  Eq.  (76),  the  factor 


/  0.2 

exp  —  sec#  J  Ka(z')dz' 


is  the  attenuation  for  the  thermal  emission  from  the  foam-covered  ocean  surface. 
For  the  US  standard  atmosphere  at  /  =  19.35  GHz,  d2  =  30  km  and  d\  =  0,  we  can 
calculate  the  attenuation  using  the  numerical  Ka  from  the  MPM.  Fig  (7)  shows  the 
attenuation  in  terms  of  the  polar  angle. 
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Figure  6:  The  brightness  temperature  of  down-going  and  up-going  waves  and  their 
comparison  with  the  approximation  formula  for  /  =  19.35  GHz. 
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Attenuation 


1 


Attenuation  of  the  Standard  Atmosphere 


Figure  7:  The  attenuation  of  the  standard  atmosphere  for  /  =  19.35  GHz.  d<i 
30  km  and  d\  =  0. 
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After  the  total  brightness  temperature  on  the  ocean  surface  is  calculated  by  [also 
from  Eq.  (32)] 

Tu(9,  <t>,  di)  =  F-  T*F(6,  <f>,  di)  +  (1  -  F)  •  T^S(9,  </>,  0),  (90) 

we  can  calculate  the  brightness  temperature  at  the  height  of  the  radiometer.  Be¬ 
tween  the  ocean  surface  and  the  radiometer,  the  atmospheric  emission  and  the  at¬ 
tenuation  must  be  taken  into  account.  At  the  radiometer  height  2  =  d2.  by  dividing 
the  constant  C  =  KB/ A2  on  both  sides  of  Eq.  (76),  the  total  brightness  temperature 
can  be  written  as 

d2  ^2 

^2  -sec $f  Ka{z”)dzn  _  -sec6  f  Ka(zf)dzf 

T*(d,  (f>,d2)  =  sec 0  j  Ra(z')T(z')e  *'  dz'  +  Tu(0,  <*>,d\)e  dl  ■  , 

di 

(91) 

where  the  first  term  is  the  thermal  emission  of  the  atmosphere  and  the  exponential 
factor  in  the  second  term  is  the  atmospheric  attenuation  for  the  wave  traveling  from 
the  ocean  surface  to  the  radiometer.  In  order  to  study  individually  the  contributions 
from  the  plain  ocean  surface,  the  foam  and  the  atmosphere,  we  re-label  the  terms 
in  Eq.  (90)  and  Eq.  (91)  as  •  - 

d2 

-s ecd  f  Ka(z')dzf 

Ts  =  (1  —  F)  •  T^S(6.  <j),  0)e  «  (92) 

d2 

-sec 6  f  Ka{z')dz' 

TF  =  F  ■TtF{0.,(b.d{)e  "1  ,  (93) 

.  d2 

7  -sec  9[Ka(z")dz" 

Ta  =  sec9  J  Ra(z')T(z')e  *'  dz'.  (94) 

di 

In  the  above  expressions,  Ts  is  the  brightness  temperature  of  the  plain  ocean  surface 
with  the  consideration  of  reflection  of  atmospheric  emission  by  water  bubbles  and 
the  attenuation  when  the  wave  travels  from  the  ocean  surface  to  the  radiometer.  T F 
is  the  brightness  temperature  of  the  foam  layer  with  the  consideration  of  scattering 
of  atmospheric  emission  and  attenuation.  TA  is  the  brightness  temperature  due  to 
the  thermal  emission  of  the  atmosphere. 

4.3  Equivalent  Polar  Angle  of  the  Spherical  Atmospheric 
Layer 

Consider  a  radiometer  placed  at  point  A  with  zenith  height  H  as  shown  in  Fig.  (8). 
We  create  a  flat  atmospheric  layer  to  approximate  the  spherical  layer  so  that  it  is 
easier  to  apply  the  geometry  for  the  RT  theory.  On  the  flat  layer,  the  equivalent 
position  of  the  radiometer  is  at  point  C,  where  we  assume  AB  =  BC.  It  can  be 
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found  that  in  the  equivalent  fiat-layer  model  of  the  atmosphere,  the  polar  angle  9' 
is  different  from  6  in  the  spherical  model.  By  projecting  the  lines  AO  and  AB  on 
the  x-  and  2-axis,  respectively,  we  find 

AB  •  sin  6  =  (R  +  H)  sin  0o,  (95) 

R  +  AB  ■  cos  9  =  ( R  +  H )  cos  90.  (96) 

By  eliminating  the  angle  90  from  Eqs.  (95)  and  (96),  it  yields 

AB  =  Vi?2 cos2  9  +  2 RH  +  H2  -  Rcos 9.  (97) 


Considering  the  triangle  BCD,  it  can  be  shown  that 

cos  9'  =  h/ab.  (98) 

Therefore  the  modified  polar  angle  in  the  equivalent  flat  atmospheric  layer  is 

H 


9'  =  cos 


-l 


VR2  cos2  9  +  2 RH  +  IP  -  R  cos  9 


(99) 


For  example,  let  the  radiometer  height  H  —  30,  000  m,  and  the  earth  radius  R  = 
6.37  x  106m.  For  the  zero  grazing  angle  (9  =  90°),  the  modified  polar  angle  can  be 
calculated  as 


9'  =  cos 


-l 


H 


s/2RH  +  H 2 


=  cos 


n/2  x  637  x  3  +  32 


=  87.2°. 


(100) 


5  Two-Scale  Model  of  Rough  Ocean  Surface 

In  the  previous  section,  we  only  considered  the  rough  ocean  surface  with  wavenumber 
less  than  5  times  the  electromagnetic  wavelength  by  defining  the  cutoff  wavenumber 
kj  —  k0/b,  where  k0  is  the  electromagnetic  wavenumber.  For  numerical  purposes, 
the  variable  l  is  defined  as  l  =  l/k  so  that  the  integration  for  l  is  from  0  to  1/kj 
instead  of  from  0  to  oo  for  k.  Therefore  numerically  kd  cannot  be  zero,  and  the 
waves  are  separated  into  large-scale  ( k  <  kd)  and  small-scale  ( k  >  kd)  categories. 
If  we  Only  consider  the  small-scale  waves,  the  total  brightness  temperature  is  not 
sensitive  to  the  value  of  kd  (kd  around  k0/ 5)  that  will  be  demonstrated  in  Fig.  (14) 
in  the  numerical  simulation  section.  The  large-scale  wave  operates  as  tilted  facet. 
The  two-scale  model  takes  care  of  the  tilted  polar  angle  on  the  local  facet  due 
to  the  large  scale  or  long  waves  of  the  ocean  surface.  In  the  two-scale  model, 
the  brightness  temperature  vector  is  calculated  by  averaging  the  local  values  over 
the  slope  distribution  of  the  large  scale  waves  [25,  18].  The  averaged  brightness 
temperature  at  the  ocean  surface  is  written  as 

OO  OO 

T  (9,  <j>,  0)  =  f  dSx  [  dSyTl(9',<j>)Ps(S'x,S'y,9),  (101) 

— oo  —  cot  6 
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Figure  8:  Geometry  of  the  spherical  atmospheric  layer. 
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Figure  9:  The  local  polar  angle  with  respect  to  the  global  looking  angle  and  the 
slope  of  the  facet. 

where  T)  {O',  <p)  is  the  brightness  temperature  of  the  local  facet  at  the  local  looking 
angle  O'.  Ps  is  the  slope  distribution  of  the  large-scale  waves  as  viewed  at  the  local 
looking  angle  O',  Sx  and  Sy  are  the  surface  slopes  along  the  global  x  and  y  axis, 
respectively,  while  S'x  and  S'y  are  the  surface  slopes  with  respect  to  the  radiometer 
observation  direction  {0',<j>).  Notice  that  Sx  is  limited  to  -  cot  0  due  to  shadowing 
by  large-scale  waves  [18].  The  transformation  of  the  global  slopes  to  the  slopes  with 
respect  to  the  radiometer  observation  direction  is  given  by 

S'x  =  Sx  cos  <*)  +  Sy  sin  6,  (102) 

S'y  =  —Sx  sirup  +  Sy  cos  <fi.  (103) 

To  calculate  the  polar  angle  O'  with  respect  to  the  local  facet  II  for  a  given 
looking  angle  0  with  respect  to  the  global  coordinate  {x.y,  z)  as  shown  in  Fig.  (9), 
we  write  the  unit  normal  vector  of  the  local  facet  as 


z  Sxx  Syy 

P+W+M’ 


and  the  unit  wave  vector  with  the  orientation  angles  0  and  p  as 

A 

k  =  x  sin  0  cos  $  +  y  sin  0  sin  4>  +  z  cos  0. 

Thus  the  polar  angle  with  respect  to  the  local  facet  is 

..  i  — Sx  sin  0 cos  (p  Sy  sin  0  sin  (p  -(-  cos  0 

0  =  cos  k  •  n  =  cos  - ,  — - . 

\/i  +  si  +  s* 


(104) 


(105) 


(106) 


The  slope  distribution  Ps  was  studied  by  Cox  and  Munk  [26]  by  measuring  the 
ocean  surface,  and  it  can  be  written  as 

p»(si,s;,f>)  =  (l  +  s1tanfl)p(s;,5;),  (107) 

26 


/ 


JL 


where  P{S'X,S')  is  assumed  to  be  a  Gaussian  distribution  function 
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2  S2  +  2SI 


with  the  function  F(S'X,  S')  defined  as 


G40 


24  l  5C4 


V  n  C22 

f +  3  +X 


9' 2 

-3L.  -  1 
52 


24  U  *2  T 


(108) 


S? 

s2 


(109) 


In  the  above  expression,  the  coefficients  are  C40  =  0.4,  C22  =  0.12,  C04  =  0.23, 
C2i  =  0.01  —  0.0086C4,,  and  C03  =  0.04  -  0.03314,,  where  Uw  is  the  wind  speed  in 
m/s. 

.  The  upwind  and  crosswind  slope  variances  are  calculated  as  follows:  .  ...  . 


K<£ 

Sl  =  J  dk  j  d<j>k 3  cos2  <f>W  (k,  <j>) , 


(110) 


Kd  *7i 

=  j  dk  J  dcpk3  sin2  <j>W  (k,  <j>) . 


(Ill) 


I11  the  numerical  evaluation  of  Eq.  (101),  the  integration  limits  of  Sx  and  Sy  are 
truncated  as  5  Su  and  5  Sc,  respectively.  The  local  brightness  temperature  T;  is 
assigned  to  be  the  one  of  the  plain  ocean  surface  Tu  as  in  Eq.  (60)  or  of  the  foam 
as  in  Eq.  (30)  based  on  the  slope  of  the  surface  S  =  \JS^  +  Sjj.  To  the  best  of 
our  knowledge,  there  is  presently  no  literature  on  the  study  of  the  foam  assignment 
according  to  the  slope  of  the  ocean  surface.  In  this  thesis,  we  assume  a  threshold 
Std  and  assign  the  foam  brightness  temperature  Tu  to  the  local  Ti  if  the  slope 
yjsl  +  S'2  >  Std •  The  slope  threshold  Std  can  be  found  by  calculating  the  integral 

within  the  area  [as  shown  in  Fig.  (10)]  of  j(Sx,  Sy)  :  S'2  +  52  <  52d,  Sx  >  —  cot#}, 
so  that  the  integral  value  is  equal  to  the  foam  coverage  F  in  Eq.  (33),  i.e. 


F  =  J  dSy  JdSxPs(S'x,S'x) 


(112) 
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Sx  =  -  cot  9 


Figure  10:  The  integration  area  for  determining  the  slope  threshold. 

6  Numerical  Results 

In  this  section,  we  calculate  the  brightness  temperature  of  the  wind-driven  ocean 
surface  and  compare  the  simulation  results  with  the  data  from  JPL’s  WINDRAD 
experiment  [14].  In  the  experiment,  a  K-band  (19.35  GHz)  radiometer  was  mounted 
on  the  NASA  DC-8  aircraft  flying  in  a  circle  at  the  height  of  30,000  ft  (9,144  nr).  The 
data  was  collected  in  November  1993  near  the  northern  Californian  coast.  During 
the  experiment,  weather  was  clear  and  there  was  a  wind  speed  of  12  m/s.  The 
brightness  temperatures  were  measured  for  the  Stokes  parameters  Tv,  Th ,  and  U 
with  the  polar  angles  of  30,  40,  and  50  degrees.  In  the  simulation,  the  frequency  of 
the  brightness  temperatures  is  /  =  19.35  MHz,  the  ocean  wind  speed  is  Uw  =  12  nr/s 
at  height  h  =  19.5  m.  The  physical  temperature  of  the  sea  water  is  T0  =  12  °C,  and 
the  salinity  is  S  =  3.5%. 

Fig.  (11),  (12)  and  (13)  show  the  simulation  results  of  the  brightness  tempera¬ 
tures  by  varying  the  azimuth  angle  4>  of  the  observation  with  nadir  looking  angle 
6  =  30°,  40°  and  50°,  respectively.  The  ocean  surface  spectrum  is  the  empirical 
formula  proposed  by  Durden  and  Vesecky  [17].  For  the  foam  layer,  the  internal 
radius  of  water  bubble  is  R  =  4.3  mm,  the  bubble  film  thickness  is  d  —  0.13  mm  and 
the  fractional  volume  of  water  bubbles  is  fv  =  0.01.  The  extinction  and  absorption 
coefficients  are  calculated  as  Ke  =  8.298  m-1  and  Ka  —  2.273  m-1.  The  permittivity 
of  the  bubble  film  is  the  same  as  of  the  ocean  water,  and  the  permittivity  of  the 
background  is  e0.  The  elevation  height  of  the  radiometer  is  30,000  ft  (9,144  m).  In 
these  figures,  the  open  circles  are  the  WINDRAD  experimental  data.  In  the  plots  of 
the  first  and  second  Stokes  parameters  Tv  and  Th,  the  lines  (from  bottom  to  top)  are 
the  plots  for  the  numerical  simulation  considering  (1)  only  the  plain  ocean  surface 
[Ts  in  Eq.  (92)],  (2)  plain  ocean  surface  plus  the  foam  emission  [TS  +  TF  in  Eq.  (92) 
and  (93)],  and  (3)  plain  ocean  surface  plus  the  foam  and  the  atmospheric  emissions 
[Ts  +  TF  +  Ta  in  Eq.  (92),  (93)  and  (94)].  In  the  plots  for  the  thiM  and  fourth 
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Figure  11:  The  brightness  temperature  of  wind-driven  ocean  surface  for  nadir  look¬ 
ing  angle  6  =  30°  and  the  lower  cutoff  wavenumber  kd  =  80  m-1. 

Stokes  parameters  U  and  V,  the  total  emissions  are  considered.  In  comparison  with 
experiment  data,  we  notice  that  both  the  foam  and  the  atmospheric  emission  are 
significant  to  correct  the  emission  of  the  plain  ocean  surface  for  Tv  and  Th- 

The  Durden- Vesecky  spectrum  is  similar  to  the  power  law  W(k)  =  a0/k4  that 
describes  the  relative  portion  of  the  large  scale  and  small  scale  roughness  of  the 
ocean  surface  by  specifying  the  lower  cutoff  wavenumber  kd.  The  smaller  kd  is.  the 
higher  the  long  waves.  However,  in  the  simulation  for  the  thermal  emission  from 
the  ocean  surface,  the  brightness  temperatures  are  obtained  by  integrating  over 
the  entire  reflected  waves  due  to  Kirchhoff’s  law.  Therefore  there  is  no  significant 
difference  between  the  collection  of  reflected  waves  from  very  long  ocean  waves  or 
from  a  flat  surface.  In  SPM,  the  zeroth  and  the  second  order  reflected  waves  are 
specular,  thus  they  include  the  dominant  reflections  from  the  long  waves.  The 
bistatic  pattern  of  the  field  scattered  by  rough  ocean  surface  is  dominated  by  small 
scale  roughness  (Bragg  scattering)  which  is  included  in  the  first  order  solution  in 
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Figure  12:  The  brightness  temperature  of  wind-driven  ocean  surface  for  nadir  look¬ 
ing  angle  6  =  40°  and  the  lower  cutoff  wavenumber  kd  —  80  m-1. 


Brightness  Temp.  [K]  _V  Brightness  Temp.  [K] 


Azimuth  Angle  [deg]  Azimuth  Angle  [deg] 


Figure  14:  The  brightness  temperature  of  wind-driven  ocean  surface  for  nadir  look¬ 
ing  angle  9  =  40°  and  the  lower  cutoff  wavenumber  kj  =  120  m-1. 

SPM.  Therefore  the  value  of  the  lower  cutoff  wavenumber  kd  is  not  sensitive  to  the 
calculation  of  the  emissivity  from  the  rough  ocean  surface.  This  can  be  demonstrated 
by  re-calculating  the  brightness  temperatures  shown  in  Figs.  (11),  (12)  and  (13)  with 
lower  cutoff  wavenumbers.  In  Figs.  (11),  (12)  and  (13),  the  lower  cutoff  wavenumber 
is  kd  =  80  m-1  which  is  about  5A  of  the  electromagnetic  wave.  Fig.  (14)  shows  the 
brightness  temperatures  for  the  polar  angle  6  =  40°  and  the  lower  cutoff  wavenumber 
kd  =  120  m-1.  Not  much  change  is  observed  for  the  result  in  comparison  with 
Fig.  (12)  for  kd  =  80  m-1. 

With  the  same  simulation  conditions  as  in  Fig.  (11)— (13),  we  compare  the  nu¬ 
merical  results  obtained  by  using  the  one-scale  and  two-scale  models  as  shown  in 
Fig.  (15)-(17).  Notice  that  the  Tv  and  Th  terms  in  these  figures  match  better  with 
the  measurement  data  in  the  two-scale  model  for  0  =  30°  and  40°,  but  there  is  an 
irregular  offset  in  Fig.  (17),  which  may  be  due  to  the  shadowing  effect  at  large  polar 
angles. 
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Figure  17:  The  comparison  of  the  one-scale  and  two-scale  models  at  6  =  50' 


7  Conclusions 


In  this  chapter,  one-scale  and  two-scale  electromagnetic  models  to  calculate  the 
brightness  temperature  of  wind-driven  ocean  with  foam  coverage  have  been  pre¬ 
sented.  The  one-scale  emissivity  model  is  the  local  thermal  contribution  by  small 
roughness  of  the  foam-covered  ocean  surface,  while  the  two-scale  model  is  the  aver¬ 
age  of  one-scale  emissivity  over  large-scale  slope  of  the  rough  ocean  surface  with  the 
weight  described  by  the  slope  distribution  function.  The  overall  brightness  tempera¬ 
ture  is  contributed  by  three  portions  —  plain  ocean  surface,  foam,  and  atmospheric 
layer.  The  interactions  between  the  different  regions  are  described  by  boundary 
conditions.  For  the  plain  ocean  surface,  the  emissivity  has  been  calculated  using 
KirchhofFs  law  by  calculating  the  reflectivity  of  rough  ocean  surface.  The  coupling 
to  the  wind  speed  is  through  the  power  spectrum  of  ocean  surface  by  which  the 
ocean  rough  surface  is  described.  The  radiative  transfer  theory  (RT)  is  used  to 
model  the  foam  layer  in  which  we  assume  the  scatterers  are  spherical  thin-film  wa¬ 
ter  bubbles.  The  closed  form  solution  of  the  RT  equations  for  the  foam  layer  is 
obtained  using  the  iterative  approach  up  to  the  first  order.  In  the  RT  equations 
for  the  atmospheric  layer,  the  scattering  due  to  the  water  vapor  and  other  gases  in 
the  air  is  ignored,  hence  a  closed  form  of  the  RT  solution  is  obtained.  The  sim¬ 
ulation  results,  both  from  one-scale  and  two-scale  models,  are  compared  with  the 
WINDRAD  experimental  data  with  good  agreements.  The  results  show  that  both 
the  one-scale  and  two-scale  models  agree  well  with  the  WINDRAD  data.  However 
the  one-scale  model  is  much  faster  in  computation  than  the  two-scale  model,  since 
few  integrals  are  involved  in  the  one- scale  model. 

It  has  to  be  pointed  out  that,  although  the  foam  model  with  water  bubbles  is 
more  realistic  than  water  particles,  more  studies  need  to  be  conducted  to  model 
the  foam  layer  more  accurately.  For  example,  it  is  suggested  the  multiple  scattering 
among  water  bubbles  to  be  considered  by  using  the  dense  medium  radiative  transfer 
theory  [15],  and  the  top  foam  surface  to  be  rough  instead  of  flat  as  considered  in 
this  thesis  for  simplification. 
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APPENDIX 


A  Phase  Matrix  for  Water  Bubbles 

Consider  an  electromagnetic  plane  wave  propagating  in  +z'  direction  and  incident 
upon  a  thin-film  water  bubble  in  (x1,  y\  z')  coordinates  with  inner  radius  R  and  film 
thickness  d  as  shown  in  Fig.  (18).  Let  the  wavenumber  in  the  core  and  outside  be  k. 
the  wavenumber  in  the  film  be  k\,  and  the  polarization  angle  be  /?.  On  y'z'  plane, 


Figure  18:  EM  scattering  by  a  bubble  in  (x'.  y'.z')  coordinates. 


in  which  m  =  ki/k,  x  =  —  Rsm(6'/2),  6  =  2nd/ A,  a  =  2nR/X,  j0(x)  and  j x(x)  are 
the  zeroth  and  first  order  spherical  Bessel’s  functions,  respectively.  The  geometry  of 
the  wave  scattering  in  {x\y\z’)  coordinates  is  general  because  of  the  symmetrical 
property  of  the  spherical  water  bubble. 
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Figure  19:  The  transformation  of  coordinates. 


In  the  (a:,  y.  z )  coordinate  system,  the  incident  wave  is  written  as  : 

Ei  =  etE0efk>f,  (A-5) 

where  e,  =  fa  or  hi,  Ee%  =  EVi  or  Ehi,  and 

ki  =  kk  =  k  ( x  sin  8i  cos  fa  +  y  sin  9i  sin  fa  +  z  cos  $i) , 

and  the  incident  angles  Oi  and  fa  are  with  respect  to  the  ( x,y,z )  coordinates  as 
shown  in  Fig.  (19).  The  polarization  vectors  are  defined  as 


fa  =  fa  x  ki  —  x  cos  9i  cos  fa  +  y  cos  6i  sin  fa  —  z  sin  6t 

hi  =  —  X  =  —x  sin  fa+y  cos  fa. 
z  x  ki 


(A-6) 


In  the  (x.  y,  z)  coordinate  system,  we  set  up  the  coordinates  (x',  y’ ,  z1)  so  that 

-  A  A 

_  ki  x  ks 

I  k  x  L 


kiXkiX  ics 

h  1 

■  h) 

1  -  k 

ki  x  hi  x  ks 

k  1 

[h  ■  k) 

1  -k 

(A-7) 


z'  =  ku 


where 


ks  —  x  sin  6S  cos  fa  +  y  sin  8S  sin  fa  +  z  cos  8S 


is  the  wave  vector  of  the  scattered  field.  Therefore  the  transformation  relation  of 
the  two  coordinate  systems  is 


or 


where 


X'  =  xan  +  yai2  +  za 13 


y'  =  xa2i  +  ya22  +  za23 
{  z'  =  xa3 1  +  ya32  +  za33 


x' 


an  = 


ai2 


®13  = 


ro 

=  A 

\ - 

N> 

1 _ 

022 

O23 

«32 

033 

031 

O33 

021 

023 

021 

031 

O22 

O32 

(A-8) 


(A-9) 


(A-10) 


a2 1  = 


121 


\/  All  +  A22  +  A| 

A22 


022  = 


23 


a23 


\/  ^21  +  A22  +  A23 
_ A23 _ 

\/  A21  +  A22  +  A23 


f  031  =  sin  cos  fa 
a32  =  sin  9i  sin  fa 
[  a33  =  cos  6i 


f  A  =  sin  6i  sin  6S  cos  (fa  —  fis)  +  cos  9X  cos  9S 
A21  =  A  sin  9i  cos  fa  —  sin  9 s  cos  &s 
A22  =  A  sin  9i  sin  fa  —  sin  9S  sin  fa 
l  A23  =  A  cos  0,- -  cos 


(All) 


(A-12) 


(A-13) 


By  using  the  identity  A  =  A  ,  the  transformation  of  a  vector  from  one  coordinate 
system  to  another  can  the  expressed  as 


E'  =  E't 


X 


=  E,tA 


=  E. 


(A-14) 


Thus  E  —  A  E'  and  E'  =  AE. 

By  writing  the  incident  field  in  (x,  y,  z)  coordinates  as 


Ei  =  eiEjh'f  =  E0  [(fa  •  x)  x  +  (fa  -y)  y  +  (fa  ■  z)  z)  e^f, 


(A- 15) 
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thus  the  transformation  of  the  incident  field  is 


'Ey 

‘(ii-x)' 

an 

<*12 

<*13 

■  (e*  -  x) ' 

El> 

=  A 

&  •  y ) 

E0  = 

a2i 

a22 

<*23 

(ei  •  y ) 

[ei,  J 

[(ei-z)\ 

.<*31 

<*32 

<*33  _ 

.  (e*  •  2) . 

(A-16) 


we  can  calculate  the  polarization  angle  P  and  scattering  angle  O'  in  the  (x',  y\  zr) 
coordinates  as  follows:  _ 

„-l  b2 


E1-  •  y' 

p  =  COS-1  r~7i  =  cos 


e: 


C' 


O'  —  cos  1  ks-  z'  =  cos  1  [sin  Oi  sin  0S  cos  (<j)$  —  (pi)  +  cos  Oi  cos  0S] , 


(A-17) 

(A-18) 


where 


bi  =  an  (e*  •  x )  +  a12  (it  ■  y)  +  ai3  (ei  •  z) ,  (A-19) 

b2  =  a2i  (ej  •  x)  +  a22  (e<  •  y)  +  a23  •  z) ,  (A-20) 

f>3  =  031  •  x)  +  a32  (ej  •  y)  +  a33  (e^  •  z) ,  (A-21) 

C  =  sjbl  +  bl  +  bi  (A-22) 


In  the  (x',  y'.  2')  coordinate  system,  we  write  the  scattered  field  as 

~E'S  =  Ee'O'  +  Eftp'  =  Eg'  (y1  cos  O'  -  z'  sin  O')  -  Eftx' 

—Eft 

—  Eg'  COS  O'  , 

— Eg'  sin  O' 


(A-23) 


eikr 

where,  by  droping  out  the  factor  —  for  scattering  coefficient  calculation, 

T1 


Eft  =  ^.Sicos/?, 

rt 

(A-24) 

Eg-  =  -- S-2  sin  p. 

A/ 

(A-25) 

The  scattered  field  in  the  (x.  y,  z)  coordinate  system  is  thus 


T 

an  a2i  a3i 

— .EV 

£5  =  A  E's  = 

<*12  a22  a32 

Eg'  COS 

<*13  <*23  033 

—Eg'  sin  O' 

’  —E^an  +  Eg'  cos  0'a2i  —  Eg'  sin  0'a3 1 
=  —Eft  cii2  +  Eg'  cos  O' a22  —  Eg-  sin  0'a32 
_ -Eftd i3  +  Eg'  cos  0'a23  -  Eg'  sin  <9'a33 

By  writing  the  polarization  vectors  of  the  scattered  field  as 


vs 


< 


hs  x  ks 
z  x  kg 

z  x  ks 


=  x  cos  0S  cos  <ps  +  y  cos  0S  sin  ps 
=  —x  sin  (j)s  +  y  cos  <f>s, 


sin  0S 


(A-26) 


(A-27) 
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the  v  and  h-components  are  obtained  as  follows: 

El  =  Es  -  vs  =  (-Epa a  +  j Eg>  cos  0'a2 1  -  E6>  sin  0'a3i)  cos  8S  cos  </>s 
+  (-E^a  12  +  -Be'  cos  (9'a22  -  Eg’  sin  6|,a32)  cos  sin  </>s 
+  —  Egi  cos  O' d2z  +  Eg/  sin  O' 0,33)  sin  0S ,  (A-28) 

££  =  ~Es-hs  =  {E^an  ~  Ee>  cos  0'a2\  +  Eg>  sin#'a31)  sin^s 

-  {E$>ai2  -  E$i  cos  0'a22  +  Eg>  sin  0'a32)  cos  <f>s.  (A-29) 

Therefore  the  scattering  coefficients  are  fvv  =  El,  fhv  =  Esh  for  e,  =  v,  and  E0  =  1; 
fvh  =  El,  fhk  =  Esh  for  ii  =  hi  and  B0  =  1.  Define  the  matrix 

L{0s,<j)s;0i,(J)i)  = 

l/w|2  \U2  Mfwflh)  -Mfwflh) 

\fkv\2  \fhh\2  ReUhvfhh)  -Im  ifhvfhh) 

2Re(fvvf*hv)  2R e(fvhfth)  Re(fvvf*hh  +  fvhflv)  -lm(fvvf*hh  -  fvhf*hv)  ’ 

21m (fvvf*hv)  21m (fvhf*hh)  lm(fvvflh  +  fvhf*hv)  Mfwfb  -  Ufa)  . 

(A-30) 


the  phase  matrix  can  be  calculated  as 

P(0s,(t)s;0i,<t>i)  =  noL(6s,<f>s:,0i,<pi), 

where  n0  is  the  number  of  bubbles  per  unit  volume. 

The  extinction  coefficient  is  derived  as  following  [15,  page  147]: 

4tt  47 r 

Ke  =  no^lm{fvv(0i,4>u0u<})i)}  = -n0j?lm{iSi} 

m2  -  1  ji(x)' 


47 r 


=  n0^-Im  |(5(m2  -  l)a 

4tt  r  2t  f,  2  2rn2  +  1 
=  n0—5a  Im  <  (m  -1)1-^- 


Jo(^)  - 


777- 


j 


The  absorption  coefficient  is 

Ks  =  Ti0  J  dfls  4>s]0U  4>i)\  +  \fhv(8s,  (ps'i  Oil  <&)[  j  ) 


(A-31) 


(A-32) 


(A-33) 


which  can  be  evaluated  numerically.  Due  to  the  symmetry  property  of  the  bubble, 
ks  is  independent  on  the  angles  0t  and  <&.  The  absorption  coefficient  is  calculated 
as 

Ka  =  Ke  -  Ks.  (A-34) 
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B  The  Parameters  in  the  Millimeter- Wave  Prop¬ 
agation  Model 


The  parameters  used  in  the  MPM  [24]  are  as  follows: 

Barometric  pressure 

P  =  P  +  e,  (B-l) 

where  p  is  dry  air  pressure  and  e  is  partial  water  vapor  pressure.  The  unit  of  the 
barometric  pressure  is  kPa. 

Temperature 

T  =  300/f,  (B-2) 


where  the  unit  of  the  temperature  is  Kelvin  (K),  and  t  is  the  inverse  temperature 
parameter.  The  typical  temperature  profile  is  given  in  [23]  and  plotted  in  Fig.  (4). 
Relative  humidity 


9  834 

RH  =  —  x  100  =  41.51-  x  lO-^, 

es  ••••  '•  t5  ■  ■■■■'■  • 


(B-3) 


where  es  is  the  saturation  pressure  over  liquid  phase. 
Dry  air  and  vapor  densities 


u  =  11.612pt, 

v  =  7.217et,  (B-4) 

where  the  unit  of  the  densities  is  g/m3. 

Complex  refractivity 

N  =  N0  +  N'(f)  +  iN"(f),  (B-5) 

where  AT0.  N'(f)  and  N"(f )  are  real  and  called  the  frequency-independent  term, 
refractive  dispersion  and  absorption,  respectively.  /  is  frequency  in  gigahertz  (GHz). 

Frequency-independent  term 

N0  =  (2.588p  +  2.39e)i  +  Nv,  (B-6) 

where  Nv  =  41. Get2  is  the  contribution  from  the  rotational  spectrum  of  water  vapor. 

Dispersion  term 

na  nb 

N'(f)  =  £  (S„F')i  +  «■;  +  £  (S^FJ),  +  K  +  K,  (B-7) 

i=l  t=l 

Absorption  term 

na  rib 

N"(f)  =  £  +  n; + £  (St.Fi'),  +  K  +  K,  (B-s) 

i= 1  i= 1 


42 


where  Sa  =  aipt3ea2^~^  and  St,  =  b\et3 -5eb2^l~^  are  the  line  strength  in  kilohertz  for 
oxygen  and  water,  respectively.  F'a  and  F6'  are  the  real  parts  of  a  line  shape  function 
in  GHz-1  which  can  be  written  explicitly  as 


X  a 


—  +  sj±-  -y)—, 

Vao  \*Ca  Vao 


(B-9) 


where  a  =  a,  b,  and 


Vao  ~  /  Vqq  +  /]  /  /  J_  J_\ 

Xa  Ya  J  Vao  \Xa  YJ 


7a/ 


(B-10) 


Xa  =  (Vao-/)2+7a> 

^a  =  (^ao  +  /)2  +  7q> 

=  (v20  +  7a)/^ao, 

7a  =  «3  (pf°'8_°4  +  l.let)  . 

76  =  &3  (pt0'8  +  4.8ef)  , 

<5a  =  a5pf06, 

<56  =  0. 


In  the  above  expressions,  voa  and  a,  {i  —  i,  2,  •  •  •  ,  6)  are  oxygen  line  parameters, 
i-ob.  and  bj  ( i  —  1,2,3)  are  water  vapor  line  parameters. 

Dry  air  continuum 

K(!)  =  oo  { [l  +  (//7o)2] -  1 }  pi2,  (B-ll) 

N;C!)=  {2oo  [70  (l  +  (//7o)2)  (l  +  (//60)2)]"‘  +  appt25)  fpt 2,  (B-12) 

where  a0  =  3.07  x  10“4,  ap  =  1.40(1  -  1.2/1 510“5)  lO”10,  and  y0  =  5.6  x  10~3(p  + 
l.le)t0-8  GHz. 

Water  vapor  continuum 

K(f)  =  b0f205et2A,  (B-13) 

K'(f)  =  (bfP  +  beet3)  fet2\  (B-14) 

where  b0  =  6.47  x  10"6,  bf  =  1.40  x  10~6  and  be  =  5.41  x  10~5. 

Hydrosol  continuum 

jv;(/)  =  2.4  x  10~W,  (B-15) 

K(f)  =  4.50to/e"(l  +  i,2),  (B-16) 

where  ij  =  (2  +  and  r  =  4.17  x  10-5te(7 13tl  ns.  e'  and  e"  are  the  real  and 
imaginary  parts  of  the  dielectric  constant  of  water.  They  can  be  calculated  using 
the  following  empirical  formulas 


e'  =  4.9  + 


185  -  113/t 

1  +  (/r)2 


(B-17) 


43 


(185-  113/Q  fr 

1  +  (fr)2 


References 

[1]  S.  H.  Yueh,  R.  Kwok,  F.  K.  Li,  S.  V.  Nghiem,  and  W.  J.  Wilson,  “Polarimetric 
passive  remote  sensing  of  ocean  wind  vectors,”  Radio  Sci,  vol.  29,  no.  4,  pp. 
799-814,  July- Aug.  1994. 

[2]  S.  H.  Yueh,  S.  V.  Nghiem,  R.  Kwok,  W.  J.  Wilson,  F.  K.  Li,  J.  T.  Johnson, 
and  J.  A.  Kong,  “Polarimetric  thermal  emission  from  periodic  water  surface,” 
Radio  Sci.,  vol.  29,  no.  1,  pp.  87-96,  Jan.-Feb.  1994. 

[3]  J.  T.  Johnson,  J.  A.  Kong,  R.  T.  Shin,  D.  H.  Staelin,  K.  O’Neill,  and  A.  W. 
Lananick,  “Third  Stokes  parameter  emission  from  a  periodic  water  surface,” 
IEEE  Trans.  Geosci.  Remote  Sensing,  vol.  31,  no.  5,  pp.  1066-1080,  Sep.  1993. 

-  [4]  J.  T.  Johnson,  J.  A.  Kong,  R.  T.  Shin,  S.  H.  Yueh,  S.  V.  Nghiem,  and  R.  Kwok, 
“Polarimetric  thermal  emission  from  rough  ocean  surfaces,”  J.  Electromagnetic 
Waves  and  Applications,  vol.  8,  no.  1,  pp.  43-59,  1994. 

[5]  G.  F.  Williams,  “Microwave  radiometry  of  the  ocean  and  the  possibility  of 
marine  wind  velocity  determination  from  satellite  observations,”  J.  Geophys. 

•-  Res.,  vol.  74,  no.  18,  pp.  4591-4594,  Aug.  1969. 

[6]  G.  F.  Williams,  “Microwave  emissivity  measurements  of  bubbles  and  foam,” 
IEEE  Trans.  Geosci.  Electron.,  vol.  GE-9,  no.  4,  pp.  221-224,  Oct.  1971. 

[7]  P.  M.  Smith,  “The  emissivity  of  sea  foam  at  19  and  37  GHz,”  IEEE  Trans. 
Geosci.  Remote  Sensing ,  vol.  26,  no.  5,  pp.  541-547,  Sep.  1988. 

[8]  A.  Stogryn,  “The  emissivity  of  sea  foam  at  microwave  frequencies,”  J.  Geophys. 
Res.,  vol.  77,  no.  9,  pp.  1658-1666,  Mar.  1972. 

[9]  P.  C.  Pandey  and  R.  K.  Kakar,  “An  empirical  microwave  emissivity  model  for 
a  foam-covered  sea,”  IEEE  J.  Oceanic  Eng.,  vol.  OE-7,  no.  3,  pp.  135-140, 
July  1982. 

[10]  J.  D.  Droppleman,  “Apparent  microwave  emissivity  of  sea  foam,”  J.  Geophys. 
Res.,  vol.  75,  no.  3,  pp.  696-698,  Jan.  1970. 

[11]  P.  W.  Rosenkranz  and  D.  H.  Staelin,  “Microwave  emissivity  of  ocean  foam  and 
its  effect  on  nadiral  radiometric  measurements,”  J.  Geophys.  Res.,  vol.  77,  no. 
33,  pp.  6528-6538,  Nov.  1972. 

[12]  X.  Z.  Huang  and  Y.  Q.  Jin,  “Scattering  and  emission  from  two-scale  randomly 
rough  sea  surface  with  foam  scatterers,”  IEE  Proc.  H  (Microwaves,  Antennas 
and  Propagat.),  vol.  142,  no.  2,  pp.  109-114,  Apr.  1995. 

[13]  E.  C.  Manahan  and  G.  MacNiocaill,  Oceanic  Whitecaps  and  Their  Role  in 
Air-Sea  Exchange  Processes,  D.  Reidel  Pub.  Co.,  Boston,  1986. 


45 


[14]  S.  H.  Yueh,  W.  J.  Wilson,  F.  K.  Li,  S.  V.  Nghiem,  and  W.  B.  Ricketts,  “Polari- 
metric  measurements  of  sea  surface  brightness  temperatures  using  an  aircraft 
K-band  radiometer,”  IEEE  Trans.  Geosci.  Remote  Sensing,  vol.  33,  no.  1,  pp. 
85-92,  Jan.  1995. 

[151  L.  Tsang,  J.  A.  Kong,  and  R.  T.  Shin.  Theory  of  Microwave  Remote  Sensinq, 
Wiley,  New  York,  1985. 

[16]  R.  T.  Shin  and  J.  A.  Kong,  “Radiative  transfer  theory  for  active  remote  sens¬ 
ing  of  two-layer  random  medium,”  in  Progress  in  Electromagnetics  Research 

■  (PIER  1),  J.  A.  Kong,  Ed.  Elsevier  Science,  New  York,  1989. 

[17]  S.  P.  Durden  and  J.  F.  Vesecky,  “A  physical  radar  cross-section  model  for  a 
wind-driven  sea  with  swell,”  IEEE  J.  Oceanic  Eng.,  vol.  OE-IO,  pp.  445-457, 
Oct.  1985. 

[18]  S.  H.  Yueh,  “Modeling  of  wind  direction  signals  in  polarimetric  sea  surface 
brightness  temperature,”  IEEE  Trans.  Geosci.  Remote  Sensing,  vol.  35,  no.  6, 
pp.  1400-1418,  Nov.  1997. 

[19]  H.  J.  Liebe,  “MPM— an  atmospheric  millimeter-wave  propagation  model,”  Int 7 
J.  Infr.  Millimeter  Waves.,  vol.  10,  no.  6,  pp.  631-650,  June  1989. 

[20]  G.  E.  Hunt.  “A  review  of  computational  techniques  for  analysing  the  transfer 
of  radiation  through  a  model  cloudy  atmosphere,”  J.  Quant.  Spectrosc.  Radiat. 
Transfer.,  vol.  11,  no.  6,  pp.  655-690,  June  1971. 

[21]  G.  Brussaard  and  P.  A.  Watson,  Atmospheric  Modeling  and  Millimeter  Wave 
Propagation,  Chapman  &  Hall,  1995. 

[22]  C.  P.  Yeang,  S.  H.  Yueh,  K.  H.  Ding,  and  J.  A.  Kong,  “Atmospheric  effect  on 
microwave  polarimetric  passive  remote  sensing  of  ocean  surfaces,”  Radio  Sci., 
vol.  34,  no.  2,  pp.  521-37,  Mar.-April.  1999. 

[23]  U.S.  Standard  Atmosphere,  National  Oceanic  and  Atmospheric  Administration, 
National  Aeronautics  and  Space  Administration,  United  States  Air  Force,  1976. 

[24]  H.  J.  Liebe,  “An  updated  model  for  millimeter  wave  propagation  in  moist  air,” 
Radio  Sci,  vol.  20,  no.  5,  pp.  1069-1089,  Sep.-Oct.  1985. 

[25]  H.  L.  Chan  and  A.  K.  Fung,  “A  theory  of  sea  scatter  at  large  incident  angles,” 
J.  Geophys.  Res.,  vol.  82,  no.  24,  pp.  3439-3444,  Aug.  1977. 

[26]  C.  S.  Cox  and  W.  H.  Munk,  “Measurement  of  the  roughness  of  the  sea  surface 
from  photographs  of  the  sun  glitter,”  J.  Opt.  Soc.  Am..,  vol.  44,  no.  11,  pp. 
838-850,  1954. 

[27]  S.  R.  Aragon  and  M.  Elwenspoek,  “Mie  scattering  from  thin  spherical  bubbles,” 
J.  Chem.  Phys.,  vol.  77,  no.  7,  pp.  3406-3413,  Oct.  1982. 


46 


OSP 


ATTACHMENT  NUMBER  1 

302 

REPORTS  AND  REPORT  DISTRIBUTION 

REPORT  TYPES 

(a)  Performance  (Technical)  Report(s)  (Include  letter  report(s))  Frequency:  Annual 

(b)  Final  Technical  Report,  issued  at  completion  of  Grant. 

NOTE:  Final  Technical  Reports  must  have  a  SF-298  accompanying  them.  Unless  otherwise 

stated  in  the  grant,  complete  Block  12a.  of  the  SF-298:  “Approved  for  Public  Release;  distribution 
is  Unlimited.” 

(c)  Final  Financial  Status  Report  (SF  269) 


(d)  Final  Patent  Report  (DD  882) 


j  REPORTS  DISTRIBUTION  || 

ADDRESSEES 

REPORT 

TYPES 

NUMBER 
OF  COPIES 

Office  of  Naval  Research 

Program  Officer  Marie  C.  Colton  ONR  321 

Ballston  Centre  Tower  One 

800  North  Quincy  Street 

Arlington,  VA  22217-5660 

(a)  &  (b) 

w/(SF-298’s) 

3 

|  Administrative  Grants  Officer 

OFFICE  OF  NAVAL  RESEARCH  REGIONAL  OFFICE  BOSTON 

★ 

495  SUMMER  STREET  ROOM  1 03 

BOSTON,  MA  02210-2109 

(cj,  (d)  &  SF-298’s 
only  for  (a)  &  (b) 

1 

Director,  Naval  Research  Laboratory 

Attn:  Code  5227 

4555  Overlook  Drive 

Washington,  DC  20375-5326 

(a)  &  (b) 

w/(SF-298's) 

1 

Defense  Technical  Information  Center 

8725  John  J.  Kingman  Road 

STE  0944 

Ft.  Belvoir,  VA  22060-6218 

(a)  &  (b) 

w/(SF-298's) 

2 

|  Office  of  Naval  Research 

Attn:  ONR  00CC1 

Ballston  Centre  Tower  One 

800  North  Quincy  Street 
|  Arlington,  VA  22217-5660 

(d) 

1 

If  the  Program  Officer  directs,  the  Grantee  shall  make  additional  distribution  of  technical  reports  in 
accordance  with  a  supplemental  distribution  list  provided  by  the  Program  Officer.  The  supplemental 
distribution  list  shall  not  exceed  250  addresses. 

For  report  types  (a)  and  (b),  send  only  a  copy  of  the  transmittal  letter  to  the  Administrative  Contracting 
Officer;  do  not  send  actual  reports  to  the  Administrative  Contracting  Officer. 


N0001 4-98-1  -0643 
(Rev.  4-96) 


